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ABSTRACT 
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2 . We present high resolution H41a hydrogen recombination line observations of the 

i 1.2' (3 pc) region surrounding Sagittarius A* (Sgr A*) at 92 GHz using the Millime- 

ter Array at the Owens Valley Radio Observatory (OVRO) with an angular resolution 
of 7" X 3" and velocity resolution of 13 km s~^. New observations of H31a, H35a, 
^ ' H41a and H44a lines were obtained using the NRAO 12-m telescope, and their rel- 

■ ative line strengths are interpreted in terms of various possible emission mechanisms. 

These NRAO 12-m measurements are the most extensive recombination line survey 
of this region to date. These OVRO data also represent the highest angular resolu- 
tion and the highest sensitivity observations of Sgr A West in hydrogen recombination 
line and continuum emission at the millimeter (mm) wavelengths. Observations of 
HCO+(J = 1 — > 0) transition at 89 GHz are obtained simultaneously with a 40% 
improved angular resolution and 4-15 times improved sensitivity over the previously 
published results, and the distribution and kinematics of the dense molecular gas in 
the circumnuclear disk (CND) are mapped and compared with those of the ionized 
gas. The line brightness ratios of the hydrogen recombination lines are consistent with 
purely spontaneous emission from Tg ~ 7000 K gas with ~ 2 x 10^ cm~^ near LTE 
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condition. A virial analysis suggests that the most prominent molecular gas clumps 
in the CND have mean densities of order lO''' cm~^, sufficient to withstand the tidal 
shear in the Galactic Center region. Therefore, these clumps may survive over several 
dynamical times, and the CND may be a dynamically stable structure. We estimate a 
total gas mass of about 3 x IO^Mq for the CND. Our combined analysis of the new high 
resolution HAla and HCO''" images and our kinematic model demonstrates a widely 
spread physical and dynamical link between the molecular gas in the CND and the 
ionized gas, including along the eastern rim of the CND where a gap was previously 
suggested. 

Subject headings: Galaxy: center — galaxies: nuclei — ISM: molecules — radio contin- 
uum: ISM — radio lines: ISM 

1. Introduction 

Studying radio recombination lines (RRLs) is one of the most direct ways of probing ionized 
gas such as HII regions. Ionized gas properties such as kinematics, electron temperatures, and the 
geometrical filling factor can be directly constrained by the recombination line flux density, line 
width, free-free continuum flux density, and size of the line-emitting regions. At centimeter (cm) 
wavelengths, RRLs are primarily stimulated emission amplifying extended non-thermal synchrotron 
emission. Therefore, distinguishing emission arising from diffuse, low density gas from that from 
clumpy, high density gas is difficult. Also, pressure broadening of the line width can complicate 
the interpretation of the cm RRL observations. For instance, cm RRL emission from luminous 
starburst systems such as Arp 220 may be dominated by internally stimulated emission from dif- 
fuse ionized gas (Anantharamaiah et al. 1993; Zhao et al. 1996) while spontaneous emission 
process is important at higher frequencies (e.g. Anantharamaiah et al. 2000; Yun et al. 2003). 
Interpretation of recombination lines at millimeter (mm) wavelengths is more straightforward since 
stimulated emission and pressure broadening become negligible in most cases. Recent studies of 
mm RRLs in Galactic H II regions (Gordon & Walmsley 1990; Martin et al. 1994; Gordon 1994) 
and starburst galaxies (Seaquist et al. 1996; Puxley et al. 1997) found higher line fluxes than 
expected from the lower frequency measurements. A steep rise in line flux density with frequency, 
characteristic of optically thin spontaneous emission, is seen in nearby starburst galaxies M82 and 
NGC 253 (Seaquist et al. 1996; Puxley et al. 1997). 

The central 1 parsec region surrounding the putative central source in our Galactic Center, 
Sgr A*, is a frequent target of multi- wavelength investigations primarily because it is the nearest 
and the best laboratory for studying the active galactic nucleus (AGN) phenomenon and its fueling 
process. Currently this region appears to be relatively gas free. However, immediately surrounding 
Sgr A* is an array of ionized gas streamers ("Sgr A West"; see review by Yusef-Zadeh et al. 2000). 
These streamers have been previously studied in the 12.8 jim. [Ne II] (Lacy et al. 1980; Serabyn & 
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Lacy 1985), Pa (Scoville et al. 2003), and in cm RRL (Schwarz et al. 1989; Roberts et al. 1991; 

Roberts & Goss 1993; Roberts, Yuscf-Zadeh, & Goss 1996), and they consist of several distinct 
kinematic features. The total mass of the ionized streamers is estimated to be less than 100 Mq in 
total. Recent IR and radio observations found evidence for neutral gas within this central "cavity" 
(Genzel et al. 1985; Jackson et al. 1993; Latvakoski et al. 1999; Herrnstein & Ho 2002), 
suggesting a more substantial presence of gas in the central parsec region. Surrounding and closely 
interacting with the ionized streamers is the dense and clumpy molecular ring commonly referred 
to as the circumnuclear disk (CND). The CND is asymmetric about Sgr A*, indicating its possible 
origin related to a gravitational capture of a passing cloud or energetic disruption of a stable ring 
(Morris &: Serabyn 1996). This ring is believed to represent a circumnuclear accretion disk feeding 
the inner parsec although some have argued that the CND is not a complete ring (see Marshall et al. 
1995; Wright et al. 2001). The ionized streamers in Sgr A West are explained as cloudlets from the 
clumpy ring that after losing momentum have fallen toward the center and thereby have stretched 
and have been ionized by photo-dissociation within the inner parsec (see Sanders 1998). Sensitive 
imaging studies in cm RRL (Roberts et al. 1991) suggest high temperature (Te = 5000 — 20, 000 
K), with a possible radial temperature gradient increasing toward Sgr A*. The large proper motions 
of the stars very close to Sgr A* suggest central dark mass of 2.5-4.1 xlO^ Mq (Genzel & Eckart 
1999; Ghez et al. 2000, 2003). 

Motivated by the detection of extragalactic mm RRLs and recent new observations of the 
Galactic Center region, we have imaged the H41q: line at high angular resolution toward Sgr A* in 
order to study the nature of the ionized streamers and their relation to the dense neutral gas in and 
around the CND. At high frequencies, most of the RRL and continuum emission comes from the 
densest ionized gas (see Brown et al. 1978; Mohan et al. 2001). One of the primary goals of this 
study is to examine the existence of ionized gas denser than that traced by cm-wavelength RRLs 
in this region. The nature of the radiation mechanism is investigated by analyzing the new H31q;, 
H35q;, H41q;, and H44a line observations obtained at the NRAO 12-m telescope at Kitt Peak, 
Arizona. We further examine whether any dense neutral gas is directly associated with the ionized 
streamers using our sensitive, high resolution HCO^ (1^0) images. Our HCO"*" (1-0) data have 
superior sensitivity and angular resolution over the previously published observations of similar 
molecular species (Glisten et al. 1987; Marr et al. 1993; Wright et al. 2001). In addition, a 
detailed analysis of the spatial distribution and kinematics of the dense molecular gas and ionized 
gas in the CND is offered using a simple kinematical model. 

2. Observations 

2.1. OVRO Millimeter Telescope Array 

The H41q! and HCO"*" (1-0) observations were made using the low, equatorial and high res- 
olution configurations of the millimeter telescope array at Owens Valley Radio Observatory. The 
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observational parameters from the multi-configuration observations conducted during the spring 
and fall 1998 arc given in Table 1. The OVRO array consists of six 10.4 meter telescopes equipped 
with cooled SIS receivers. Interferometric measurements are made using a total of 45 independent 
baselines in order to maximize the sampling of the uv-plane and to achieve the highest imaging 
dynamic range possible. In the high resolution configuration the longest baseline used was 220 
meters. The spectrometer modules were arranged so that each line was covered with a total of 60 
channels at 4 MHz (13 km s"^) resolution centered at the velocity of Vlsb. = 0.0 km s~^. By 
setting the LO frequency at 90.5 GHz, H41a line at 92 GHz (upper sideband) and HCO"^ (1^0) 
line at 89 GHz (lower sideband) could be observed simultaneously. The primary beam (field of 
view) is about 85" at these frequencies. A nearby calibrator NRAO 530 was observed every 15 
minutes to track the atmosphere and instrumental gain. The amplitude scaling is derived from the 
observations of 3C 273, whose flux density is monitored daily relative to Uranus and Neptune with 
an accuracy of about 20%. 

The data were reduced using the OVRO data reduction program MMA (Scoville et al. 1993). 

Sgr A* is a low declination source, S ~ —29°, and the synthesized beam is highly elongated in the N- 
S direction. The resulting beam size is 6.95" x 3.47" (FWHM). The data were imaged and analyzed 
using DIFMAP (Pearson et al. 1994) and Astronomical Image Processing System (AIPS) provided 
by the National Radio Astronomy Observatory (NRAO). Of the total 60 channels in the spectra, 
15 line-free channels (60 MHz total) were used to construct the continuum image. Unresolved 
strong continuum source Sgr A* (2.0 ± 0.3 Jy beam~^) and the ionized streamers were used for 
self-calibration in DIFMAP, and the continuum was subtracted from the spectral line channels. All 
spectral and continuum images were corrected for primary beam attenuation. Subsequent analysis 
of the spectral line data cubes were conducted using AIPS and Groningen Image Processing System 
(GIPSY). 



2.2. NRAO 12-m Telescope 

We have obtained the H31a (210.502 GHz), H35a (147.047 GHz), H41a (92.034 GHz), and 
H44q; (74.645 GHz) hydrogen recombination line spectra of Sgr A West using the NRAO 12-m 
telescope at Kitt Peak during May 1998. The NRAO 12-m telescope is equipped with single 
sideband SIS receivers, and typical system temperature ranged between 250 and 400 K (SSB), 
despite the low elevation of the source. The beam size is inversely proportional to the frequency 
and is 69" at 92 GHz. Both the filter banks (500 MHz at 2 MHz resolution) and hybrid digital 
spectrometer (600 MHz at 0.78 MHz resolution) were used to record the data. A wide range of 
OFF positions were examined, mostly along the minima in CO and HON surveys of the Galactic 
Center region. After an extensive investigation, it was determined that the hydrogen recombination 
line emission is highly localized around Sgr A*, and the final spectra presented here were obtained 
with a simple azimuth offset of ±3'. 

From the observatory report on the telescope efficiency measurements, we adopt the corrected 
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main beam efficiency r/^ of 0.48, 0.85, 0.90, and 0.95 for the four observed frequencies, respectively 
(Mangum 2000). These efficiencies are used to convert the observed brightness temperature to 
the main beam temperature T^b and then to flux density. The telescope efficiency measurements 
and the subsequent corrections are inherently uncertain, and we adopt an absolute uncertainty of 
20% for the measured flux density. 



The 92 GHz continuum image of the Sgr A complex at the resolution of 6.95" x 3.47" (PA = 
—5°) is shown in Figure 1. The major components of the Sgr A complex detected are the spiral 
structure known as Sgr A West and the compact central source Sgr A*, marked by a cross. The 
bright central source, Sgr A*, has peak flux density of 1.95 ± 0.30 Jy beam~^. The total 92 GHz 
continuum flux density detected inside the 85" diameter region is ~12 Jy. Some of the extended 
continuum emission may be missed by our observations as the 19-pointing mosaic image by Wright 
et al. (2001) reported a total flux of ~15 Jy covering a slightly larger area while Wright & Backer 
(1993) reported a total flux of 27 Jy within a 2' radius of Sgr A* from a combination of BIMA and 
Nobeyama 45-m telescope data. 

The major continuum features of Sgr A West are identified in Figure 1, and their observed and 
derived properties are tabulated in Table 2. The east-west oriented filamentary structure located 
immediately south of Sgr A* is the bar of Sgr A West. The total continuum flux at 92 GHz from 
the bar is about 0.9 Jy with a peak brightness of 0.7 Jy beam~^. The bar appears to connect 
to the western arc and the eastern arc on each ends. Running almost perpendicular to the bar is 
the northern arm with higher total integrated flux density of about 1.1 Jy. The continuum flux is 
higher toward the southern part of the northern arm. To the west of the bar is the western arc that 
curves south. It also extends north almost reaching the northern arm. The continuum flux from 
the western arc is relatively weak. The eastern end of the bar curves along the eastern arm but is 
not detected further north. The morphology of the ionized gas around the Sgr A complex shown 
in Figure 1 is remarkably similar to the 8.3 GHz continuum image by Roberts &; Goss (1993). 

The peak flux density measured outside Sgr A* range between 0.1 — 0.75 Jy beam~^ (see Ta- 
ble 2), which corresponds to peak brightness temperature of 1.1 — 6.0 K. Assuming this continuum 
is largely thermal free- free emission from 7000 K gas (see § 4.2), one can infer an optical depth of 
Tc < lO^"' if the emitting region fills the beam with a unit filling factor. Even if the filling factor 
is not near unity, this inference of optically thin emission should be secure unless the continuum 
emission is entirely concentrated into discrete clumps with a diameter of 0.005 pc (1000 AU) or 
smaller. Continuum optical depth is given by 



3. 



RESULTS 



3.1. 



Continuum at 92 GHz 



Tc = 0.0824 EM T-^-^^Vq]j\ 



(1) 
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Assuming Tg = 7000 K, the emission measure {EM = J nidi) required to produce continuum 
opacity of 10^"^ at 92 GHz is 2.5 x lO''' cm~^ pc. If the thickness of the filaments along the line of 
sight is comparable to the apparent width [l ~ 0.1 pc), then the inferred mean electron density is 
rie ~ 2 X 10^ cm~^ (A^g ~ 5 x 10^^ cm~^). If the ionized gas is clumped at much smaller scales, then 
the electron density would be correspondingly larger. In comparison, Lutz et al. (1996) derived 
rie ~ 3000 cm~^ from their analysis of fine-structure line ratios. The consideration of collisional 
de-excitation of the observed [Ne II] emission places an upper limit on the electron density of about 
10^ cm~^ (Jackson et al. 1993). Therefore the geometry of the ionized gas must be fairly smooth 
over 0.01-0.1 pc scales. 

3.2. H41a 

The velocity integrated line flux density ("moment zero") image, shown in Figure 2, was 
generated by summing the continuum subtracted H41q; image cube over all velocity channels with 
line emission (from —293 km s~^ to -1-221 km s~^; see Figure 3). The cross marks the position of 
Sgr A*. The main features seen in the line emission map are very similar to the 92 GHz continuum 

map shown in Figure 1, except for the absence of the bright, compact source corresponding to Sgr 
A*. The emission from the western arc is weaker than seen in continuum. The northern arm is the 
brightest continuum feature besides Sgr A*, and the brightest H41a emission region also occurs in 
the northern arm just northeast of Sgr A*. The central bar shows a wide range of emission intensity 
and structures. 

A summary of the measured parameters for all major features are given in Table 3. The sizes 
tabulated for the listed regions have been calculated assuming 1' = 2.5 pc (Reid 1989). Unlike the 
92 GHz continuum that peaks at Sgr A*, the brightest feature in H41a occurs in the northern arm 
with a peak brightness temperature of ATpeafc = 9.4 it 2.2 K. The western arc, however, appears 
relatively weak in H41a line. The peak line brightness for all H41a are < 10 K, and thus the H41a 
emission is also optically thin, tl ^ 10~^ (see below). 

The spatially averaged H41a line profiles for the six prominent emission regions (identified 
with boxes) are shown in Figure 4. The GIPSY tasks PROFIL and PROFIT were used to calculate 
the line-to-continuum (L/C) ratios and to fit a Gaussian profile to the data using a least squares 
procedure. The crosses indicate the data points, the solid lines are the Gaussian profiles, and 
the dashed lines represent the residuals. Calculations of Gaussian profiles yield four parameters: 
peak line-to-continuum flux density ratio, velocity integrated line intensity, central velocity, and 
the line width. These parameters are summarized in Table 4 along with the derived LTE electron 
temperature derived from the L/C ratios (see further discussions in § 4.2). A typical L/C ratio 
along the northern arm is around 0.4. The L/C ratio is the largest (~ 0.9) in the western arc while 
the L/C ratio in the bar drops from 0.6 to 0.03 within the 2 to 4 beam radii of Sgr A*. 

The kinematics of the ionized gas traced by H41q; shown in Figures 3 and 4 are generally 
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consistent with that of rotation about Sgr A*, similar to that of the CND, in broad agreement with 
previous studies of SgrA West in RRLs (Roberts et al. 1991; Roberts & Goss 1993) and in 12.8 
/im [Nc II] (Serabyn &: Lacy 1985). The central bar spans a very wide range of velocity, between 
+33 km and —293 km s~^ (see Figure 3), much larger than those associated with the northern 
arm (between +163 km s~^ and +7 km s~^) and the western arc (between —7 km s"^ and —111 
km s"-*^). This large velocity gradient suggests that this gas is located deeper in the gravitational 
potential, much closer to the central mass concentration, than the rest of the ionized gas. Citing 
the relatively cool dust temperature along the bar, Latvakoski et al. (1999) have argued that the 
nearness of the bar to Sgr A* is largely in projection only, but this large spread in velocity and the 
correspondingly large velocity gradient requires that this gas feature has to be located much closer 
to the central mass concentration than the northern arm and the western arc. Since the CND is 
already viewed with a large inclination [i ~ 60°), placing a gas streamer in a more edge-on orbit can 
increase the line-of-sight velocity only by about 25%. Therefore a viewing geometry alone cannot 
account for the large observed velocity gradient. Neither the angular resolution nor the overall S/N 
of the new H41a data arc better than the previous VLA H92a measurements, and the main new 
insight obtained is that the ionized gas traced in H41a emission share the same kinematics as the 
gas responsible for the H92a emission. 

The H41q; line widths of the ionized streamers range between 43 km s~^ and 347 km s^^ (see 
Table 4). The broadest line is seen toward Sgr A* with AV = 347 ± 48 km s'^ (see Figure 4, Box 
4), but the S/N of this measurement is quite small in each individual channels. The next broadest 
line occurs in the vicinity of bright H41a region in the northern arm with AV = 213 ± 18 km s^^ 
(Box 6 in Figure 4). The narrowest line is seen toward the southern end of the western arc (Box 3 
in Figure 4) with AV = 43 ± 5 km s"^ Using 2" resolution VLA H92a images, Roberts & Goss 
(1993) have shown that crowding of distinct kinematical components along the line of sight can 
make the emission lines appear much broader than the intrinsic line width of about 50 km s~^. 
Combined with the large velocity gradient expected near the central mass concentration, the large 
line width we derive from H41a; emission must be primarily the result of spatial averaging resulting 
from our large observing beam of 7.0" x 3.5". A typical line width measured in each pixel is between 
40-75 km s~^. Given the low S/N of the data, it is not physically meaningful to deconvolve the 
individual kinematic components by fitting multiple Gaussians to the observed spectra. There is 
a clear trend of increasing L/C ratio with a decreasing line width. Such a trend is expected if the 
line and continuum emission arises from the same ionized gas clouds and the emitting gas has a 
narrow range of electron temperature (see § 4.2). 

3.3. HCO+(J = 1^0) 

Figure 5 shows the high resolution image of the molecular ring in HCO"*" emission (grey-scale) 
superimposed over the H41a (contour) emission. The circumnuclear ring traced in HCO+ appears 
more complete than seen in the previous HCN (Glisten et al. 1987) and HCO"'" (Marr et al. 



-8- 



1992) images obtained using the BIMA array. More recently Wright et al. (2001) imaged HCN, 

HCO"*", and H42q; in a wider field of view using a 19-pointing mosaic, offering a much better overall 
picture of the molecular gas distribution in the 10 pc region surrounding the Sgr A*. However, the 
latest BIMA data by Wright et al. are not much better than the older BIMA data in resolution 
and sensitivity. These different images agree with each other broadly, but significant differences 
also exist in the details. One important difference is that the angular resolution of our new OVRO 
data is 40% better than those of Giistcn ct al. and Marr et al., and our sensitivity is also 4 and 15 
times improved, respectively. The BIMA array is more sensitive to faint extended emission, and its 
data reduction procedure is tuned to emphasize such features. In contrast, our moment analysis 
is more sensitive to identifying bright emission clumps, and this accounts for the more complete 
appearance of the CND in Figure 5, consisting of many discrete features. 

Determining whether any dense neutral gas is associated with the ionized streamers surrounding 
Sgr A* is one of the major objectives for our study. While the new observations have revealed clear 
evidence for a close physical link between ionized and neutral molecular gas (as discussed later in 
§4.4), no HCO"'' emission is found spatially and kinematically coincident with the ionized streamers 
traced in H41a (See Figure 3). Actually HCO^ absorption is seen in many channels, particularly 
at blueshifted velocities, but most of these features show little kinematic relation to the H41a 
emission. The spatial distribution of HCO"*" absorption resembles the entire continuum structure, 
and therefore the absorbing gas must be located in the foreground, associated with the 100 pc scale 
molecular disk imaged by Tsuboi et al. (1999) and others. Among the redshifted velocity channels 
where little foreground absorption is present (e.g., +117 to +195 km/s), there is little evidence for 
HCO^ absorption at the locations of the H41q; and continuum emission peaks. In these channels, 
the upper limit to the optical depth for HCO^ absorption is about 0.3 (3cj) at a spectral resolution 
of 13 km s~^. This limit corresponds to an upper limit on the HCO'*' column density of 4 x 10^^ 
cm~^. For an assumed HCO"^ abundance of 10~^, this translates into an H2 column density of 
4 X 10^^ cm~^. If the HCO+ abundance appropriate for a translucent cloud (where ion chemistry 
is important, see Hogerheijde et al. 1995) is adopted, the inferred upper limit on the H2 column 
density becomes as little as ~ 10^^ cm^'^. These upper limits arc close to or smaller than the column 
density of Nh ~ (0.5 — 1.0) x 10^^ cm~^ inferred from the [O I]/[C II] far-IR line ratio (Jackson 
et al. 1993) and the mid-IR color analysis (Latvakoski et al. 1999). Using the calculations by 
Sternberg &; Dalgarno (1989), Jackson et al. (1993) have proposed that most of the molecules in 
the neutral gas inside the central cavity with Nh < 10^^ cm~^ are photo-dissociated and that the 
neutral gas is primarily atomic. Our HCO"'' absorption measurements are consistent with these 
conclusions, and we conclude independently that Sgr A West filaments are either highly ionized or 
are largely atomic. 
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3.4. U31a, H35a, H41a, and H44q; 

The hydrogen recombination hne spectra obtained at the NRAO 12-m telescope are shown 
in Figure 6, and the observing frequency and the measured hne flux are summarized in Table 6. 
The mm frequency window is rich in spectral transitions, and studying faint emission lines in the 
Galactic Center region can be challenging. The H35a line is not shown because it is completely 
swamped by the much stronger CS J = 3 — >■ 2 line at 146.969 GHz, and no useful information could 
be derived. 

The H44a spectrum, shown at the bottom of Figure 6, is the least confused and the highest 
S/N measurement among all transitions observed at the 12-m telescope. The observed line profile 
closely resembles the H41a spectrum derived from the OVRO interferometer data (shown with a 
solid line in the middle of Figure 6). The total velocity integrated line flux density is 232 it 46 Jy 
km within the 84" beam of the 12-m telescope at 74.645 GHz. 

The H41q; spectrum obtained using the NRAO 12m telescope, shown with a dotted line in the 
middle of Figure 6, shows a bright, narrow feature near -1-200 km s~^. This is not seen in the H41q; 
spectrum obtained from the OVRO data (shown in solid line). We have mapped H41q; emission in 
a 2' by 3' region centered on Sgr A* in order to examine the spatial extent of this spectral feature. 
The broad emission feature seen in the OVRO data and in the H44a spectrum is spatially confined 
only to the central position while the bright, narrow feature appears everywhere with little velocity 
gradient. The narrow line width, little velocity gradient, and wide-spread distribution strongly 
indicate that this narrow component is a Galactic foreground emission. A likely identification is 
CH3CN, which has transitions at 91.959, 91.970, 91.980, 91.985, and 91.987 GHz. The total velocity 
integrated H41a line flux density from the OVRO data is 152 it 30 Jy km s~^. This is a lower limit 
since some extended flux may be missed by the interferometer (see below). 

The HSlcK spectrum shown at the top of Figure 6 is obtained by summing 8 independent 
spectra covering the whole Sgr A West complex. The telescope beam at 210.502 GHz is about 
30" and spatially resolves the Sgr A West structure. Therefore a 3 x 3 grid map was obtained at 
full beam sampling. The NE grid point was not observed, and the 8 spectra are summed together 
to derive the integrated H31a spectrum. The H41a channel maps (Figure 3) suggest that our 
H31a spectrum may be missing some redshifted emission at F = 50 - 150 km s~^. Nevertheless, 
the summed spectrum closely resembles the H41q; and H44q; spectra, and we made no attempt to 
correct for any missing flux. 

4. DISCUSSION 

4.1. Massive Compact Molecular Clumps in the CND 

The nine brightest HCO"*" clumps within the CND are identified as localized peaks of velocity 
integrated HCO"^ fiux density map (Figure 5), and their physical sizes and emission characteristics 



-10- 



are summarized in Table 5. The peak line brightness temperature, averaged over the 7" x 3.5" (0.29 
pc X 0.15 pc) beam, ranges between 10 and 30 K. In comparison, the gas temperature derived from 
NH3 is 60-70 K (Coil & Ho 2000) while the dust temperature of about 75 K is derived by Latvakoski 
et al. (1999) from the far-infrared observations. Since HCO"*" emission is strongly absorbed by the 
foreground gas along the line of sight and the beam filling factor is likely less than unity, the gas 
temperature inferred from the observed line brightness is strictly a lower limit. The measured peak 
HCO"^ line temperature is comparable to the values found by earlier BIMA studies (see § 3.3) with 
a lower angular resolution, and one can safely deduce that the optically thick molecular gas nearly 
fills the beam at ~ 0.1 pc scales in these CND clumps. 

We have computed virial masses of the individual HCO+ emitting clouds as, 

M^ir = 250 rAV^ Mq, (2) 

where r is the radius of the cloud in pc and Ay is the line width in km s^^. If these individual 
molecular gas clumps arc sclf-gravitating, then their masses range between (0.3 - 4.5) x 
(see Table 5). Assuming a spherical geometry, the mean mass density {p = Myir / ^ttv'^) is between 
(1.3 — 10) X 10"^'^ g cm~^ and mean H2 molecule density of n-vir = (0.4 - 3.0) x lO''' cm~^. These 
densities are significantly larger than the values inferred from the excitation analysis (n ~ 10^~^ 
cm~^, Genzel ct al. 1985; Marr et al. 1993; Marshall et al. 1995). For a central source with a 
mass of 4.1 ±0.6 x W^Mq (Ghez et al. 2003), the critical mass density required to survive the tidal 
shear at a radius of 1.5 pc is pudal = 4.3 x 10"^'' g cm""^ or nudal = 2.5 x 10^ cm~^. Contrary to the 
previous studies using excitation analysis, we conclude that the mean gas density inferred from the 
virial analysis is close to or exceeding the critical density for surviving the tidal shear if they are 
self-gravitating. Also, the CND may well be a structure that can survive several dynamical time 
scales (> 10^ years). The total gas mass for the CND inferred from our virial analysis is around 
3 X IO^Mq, which is only a few per cent of the total dynamical mass for the central one parsec 
region. The virial mass of each of these clumps are about two orders of magnitude larger than that 
of the ionized gas associated with the whole ionized spiral (Sgr A West). Therefore, an infall of 
even a much smaller gas clump is sufficient to account for the observed ionized streamers. The new 
higher resolution 10-ficId mosaic observations of HCO"'" and HCN emission by Christopher et al. 
(2004) suggest even higher mean density (n ^ 10^ cm~^) for these dense neutral clumps. 



4.2. Electron Temperature in Sgr A West 

The 92 GHz continuum around Sgr A is primarily free-free emission, and a good spatial 
correspondence between the continuum (Fig. 1) and the H41a line (Fig. 2) emission is naturally 
expected. Despite the similar morphology, however, a comparison suggests that the line emission 
is significantly clumpier than the continuum, especially in the western arc. The central bar is also 
clumpier in the line emission. In addition, in the northern arm the line emission appears brighter 
than anywhere else with an increasing trend toward Sgr A*. These differences may not be intrinsic 
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to the source. A lower signal-to-noise ratio in the line image probably accounts for some of the 
differences. In addition, the continuum image suffers more from the missing flux problem than the 
line images because of a larger effective source size.^ 

For optically thin gas in LTE emitting thermal continuum and negligible pressure broadening, 
LTE electron temperature T* can be derived from the line-to-continuum ratio as, 

r; = (6943."^^)«- (3) 

where u is the observing frequency in GHz, /SV is the velocity width (FWHM) in km s~^, Sc and Sl 
are continuum and line flux densities, respectively, and is the singly ionized helium abundance 
(Brown et al. 1978; Roberts & Goss 1993). In the case of Sgr A West the aforementioned 
assumptions may be reasonably made along with the interpretation that the continuum at 92 GHz 
is largely free-free emission. The singly ionized helium abundance Y'^ of 10% is assumed for Sgr A 
West (Roberts & Goss 1993). The H41q! line widths listed in Table 4 are larger than the typical 
H92q; line widths of 40-50 km s"-*^ reported by Roberts & Goss (1993), and this is the direct 
consequence of a poorer angular resolution of the H41q! observations. The resulting blending of 
emission features could conspire to produce unusual L/C ratios if gas temperature and density vary 
significantly along the ionized streamers. The high angular resolution observations by Roberts &: 
Goss (1993) have shown that this is not the case, and the poorer angular resolution of our H41q; 
observations should not affect our Tg analysis in most cases. 

An electron temperature map can be derived from the 92 GHz continuum map (Fig. 1), the 
velocity integrated H41a map (Fig. 2), and the relation given in Eq. 3. The derived electron 
temperature map, shown in Figure 7, reveals a fairly uniform distribution of Tg around 7000 K 
throughout the Sgr A West complex except near Sgr A* (see Table 4). In general T* derived at 
higher frequencies tends to be larger because of diminishing contribution from stimulated emission. 
However, electron temperature derived with H41q: line at 92 GHz agrees well with the derivation 
using H92q; line at 8.3 GHz (7000±500 K, Roberts & Goss 1993), suggesting that the contribution 
from stimulated emission is negligible even frequencies as low as 8.3 GHz. Detecting little H41q; 
emission in the direction of the bright continuum source Sgr A* is another clear evidence that 
stimulated emission is not important at 92 GHz (see § 4.3). 

The higher electron temperature derived along the central bar just south of Sgr A* has previ- 
ously been noted by Schwarz et al. (1989) and Roberts et al. (1991), and whether this is a real 
effect or an artifact of data analysis has been a subject of recent debates. Citing the discovery of 
highly blueshifted line emission near —250 km s~^ missed by these earlier studies, Roberts &; Goss 
(1993) have suggested that the unusually low L/C ratio and thus the inferences of higher electron 
temperature may be erroneous. In general, a lower than expected L/C ratio can result from either 



^ A line emitting area is significantly smaller in size than the continuum emitting region in individual channel maps 
for an object with a large velocity gradient such as a rotating disk. 
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(1) missing line flux or (2) extra continuum flux. The explanation by Roberts & Goss (1993) is 
the former, but a careful analysis using the actual measurements of the extremely blueshifted line 
by Roberts, Yusef-Zadeh, & Goss (1996) revealed that the missed flux alone cannot fully account 
for the low L/C ratio. In addition, Roberts, Yusef-Zadeh, k. Goss (1996) noted the presence of 
two optically thick continuum sources IRS 2 and IRS 13 whose additional continuum flux should 
contribute to the unusually low L/C ratio (and higher Tg) observed. A broad coincidence of this 
position with that of the "mini-cavity" (Yuscf-Zadch et al. 1990) and its possible impact on 
the calculation of T* are discussed in detail by Roberts & Goss (1993). Presence of significant 
structures at scales smaller than our synthesized beam violates the simple geometrical assumption 
going into Eq. 3 (also see Brown et al. 1978), and this may account for part of the observed local 
"enhancement" in T*. 

The local enhancement in electron temperature near the location of IRS 2 and IRS 13 is 
clearly seen in Figure 7 in the same manner as the seen in Figure 14 of Roberts & Goss (1993). 
The persistence of this feature at 92 GHz implies that the presence of the two continuum sources 
as suggested by Roberts, Yusef-Zadeh, & Goss (1996) does not offer a full explanation for the 
observed low L/C ratio. Roberts, Yuscf-Zadch, & Goss (1996) estimated an optical depth of ~ 0.8 
and ~ 1.1 at 8.3 GHz for the continuum sources IRS 2 and IRS 13 by comparing their 8.3 GHz 
image with the 5 and 15 GHz VLA continuum images. Free-free optical depth scales as r oc 
(see Eq. 1), and the continuum optical depth for IRS 2 and IRS 13 should be < 0.01 at 92 GHz. 
In comparison, optically thin free-free emission associated with the ionized filament should scale 
only as Sy oc u^^ "^. As a result, the continuum contribution from IRS 2 and IRS 13 shoTild be 
substantially reduced at 92 GHz, and the L/C ratio should refiect mostly the optically thin emission 
associated with the ionized gas filaments in SgrA West. However, the observed L/C ratio at 92 
GHz is still noticeably lower at this location. The area of a lower L/C ratio is not centered on these 
two continuum sources either. If continuum opacities for these two sources are much larger than 
estimated by Roberts et al. (i.e., optically thick even at 92 GHz), then the explanation of IRS 2 
and IRS 13 as additional sources of continuum is still plausible, although the spatial distribution is 
still a problem. In summary, the apparent enhancement in Tg (a low L/C ratio) along the central 
bar just southwest of Sgr A* is still not fully explained. 



4.3. Nature of the Recombination Line Emission 

The flux density ratio between H41a and H92a lines is shown in Figure 8. The overall structure 
of the northern arm, the bar, western arc, and eastern arm is almost identical at both frequencies, 
and the line ratio is nearly constant, SH92a/SH4ia ~ 0.15. This similarity is remarkable given the 
order of magnitude difference in the frequencies of these lines. A slightly higher ratio is seen along 
the western arc and toward the top where the western arc arches into the northern arm, suggesting 
brighter emission at centimeter wavelengths. This trend of comparatively brighter emission at 
longer wavelength was also noted in the comparison of H92a to Br7 ratios by Roberts &: Goss 
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(1993). Roberts &; Goss offered extinction of Br7 emission by dust in the CND as a possible 
explanation, but a higher Suma/ Suiia ratio would require a different explanation. Since the 
derived T* distribution is fairly constant (see Fig. 7 and Table 4), the enhanced SH92a/ Snua ratio 
is not likely caused by any changes in the physical characteristics of the gas. The comparison 
of the total integrated HAla line flux with those of other transitions measured with the NRAO 
12-m telescope suggest that some extended flux is missed by the interferometer (see § 3.4). The 
shortest baseline in the OVRO observations of 3100 A is about 4 times larger than that for the VLA 
observations (830 A), and the greater missing flux in the OVRO data offers a plausible explanation 
for the apparent enhanced SH92a/ Suiia ratio along the western arc. 

The analysis of the 92 GHz continuum (§3.1) and the H41a emission (§3.2) suggests high 
electron density (ne ~ 2 x 10^ cm^'^) and low optical depth (r < 10^'^) for the ionized streamers. 
Under these conditions, both internal and external stimulated emission should be negligible, and 
optically thin spontaneous emission should dominate the recombination line emission at millimeter 
wavelengths. For the optically thin case near LTE, Eq. 3 implies that the velocity integrated line 
flux density Sl/^V should increase linearly with frequency since Sc cc v^^'^ for thermal free-free 
emission.^ A plot of the velocity integrated line flux density for H31q!, H41q;, H44a, and H92a lines 
is shown as a function of frequency shown in Figure 9, and indeed a nearly linear trend is seen. 
Excluding the H41a measurement from the OVRO (which is a lower limit because of the missing 
flux), a best fit power-law slope is a = -1-0.95 ± 0.18, which is nearly exactly the theoretically 
expected value. 

The assumption of "near LTE" is not strictly correct, however. A general expression for 
recombination line brightness in comparison to the LTE value is 

Tl = nbn{l - ^TcPn) (4) 

where T'^ is the LTE line brightness, 6„ is the departure coefficient, is continuum optical depth, 
and Pn is the effective line absorption coefficient accounting for non-LTE stimulated emission in 
the line. For the range of physical conditions we derive for the Sgr A West (Te ~ 7000 K, rie ~ 10^ 
cm~^), the computed values of these parameters are 6„ ~ 0.7-0.9 and (1 — /3„) ~ 25 (Brocklehurst 
1970; Walmsley 1990). Although the non-LTE factor in the line absorption coefficient /?„ is large, 
it has little impact on the observed line brightness because a small optical depth (tc < 10~^, as 
shown in § 3.1) makes the second term in Eq. 4 negligible, i.e. ^Tc/9„ 0. The net effect is 
that the observed recombination line brightness is redTiccd by a factor of 6,„ from the LTE value 
(i.e., Tl/T^ = hn)- The linear slope of the the observed power-law relation between the velocity 
integrated line flux density and frequency is unaffected as long as \TcPn ^ 0- As a result, the 
average flux density of the recombination line should increase linearly with frequency even when 
some non-LTE effects are present as long as spontaneous emission dominates the radiation process. 



^Using slightly different assumptions, Seaquist et al. (1996) derive a relation showing Sl^V oc v'^^ in their 
Eq. A9. However, the factor, which they treated as a "constant" term going from their Eq. A8 to Eq. A9, actually 
behaves like {vuna oc [1/n^ — l/{n+ 1)^]), and this too leads to a linear relation, Sl^V oc z/"*"^. 
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The comparison of the H41q: and H92a spectra obtained toward Sgr A* (Figure 10) offers a 
definitive proof that the narrow H92a hne seen near V = +50 km s~^ by Roberts & Goss (1993) 
is indeed stimulated emission. Citing the nonthermal nature of the Sgr A* continuum emission and 
the low observed L/C ratio, Roberts &: Goss suggested that this H92q: emission must be stimulated 
emission arising from diffuse foreground ionized gas. The H92a spectrum (shown in dotted line) 
and the H41q; spectrum (shown in solid line) have essentially identical shape at a position along 
the Northern Arm, 7" offset to the east of Sgr A*, suggesting the same ionized gas is responsible 
for both spectral features. At the location of Sgr A*, however, H92a spectrum shows a distinct 
narrow component near V = +50 km s~^, not seen in the H41q spectrum. As shown in Eq. 1, the 
optical depth for H41a line is expected to be more than 100 times smaller than the H92q; line, and 
H41a is expected to be smaller by this large factor for optically thin stimulated emission. 



4.4. Physical Link between Neutral and Ionized Gas 

A review of the recent literature yields a coherent but complex picture for the composition, 
geometry, and kinematics of the neutral and ionized gas found in the central few parsec of our 
Galaxy. The CND, traced in HCO^ and HCN, has a ring-like morphology with an inner radius 
of about 30" (1.25 pc) and an outer radius of at least 45-60" (1.8-2.5 pc) and is physically and 
kinematically linked with a larger scale disk with an extent of 10 pc or more (Wright ct al. 2001; 
McGary et al. 2001; Christopher et al. 2004). It may be made of a multi-phase medium as the 
estimates of gas temperature (50-250 K) and density (10^~^ cm~^ and > 10^ cm~^ in dense clumps, 
see § 3.3) range widely (Marr et al. 1993; Jackson et al. 1993; Marshall et al. 1995; Latvakoski 
et al. 1999). Within it lies the "central cavity" which actually contains the ionized streamers Sgr 
A West, made of about 100 AIq of rsj IQ'* cm^^ gas at Tg ^ 7000 K and as much as 1000 Mq of 
neutral gas which is thought to be mostly atomic with nn 10^ cm""^ and T = 100-300 K (Genzel 
et al. 1985; Roberts & Goss 1993; Jackson et al. 1993; Telesco et al. 1996; Latvakoski et al. 
1999, and see § 3.3). By obtaining the most sensitive and the highest angular resolution images of 
H41a, 92 GHz continuum, and HCO+ (1-0) to date simultaneously, we are able to shed some new 
light on the physical link between the neutral and ionized gas in this region. 

The clearest demonstration of the physical connection between the neutral and ionized gas 
is the presence of recombination line peaks along the western arc, just inside of the dense gas 
clumps within the CND as shown by the overlay of the velocity integrated H41a and HCO^ images 
(see Figure 5). This radial offset between the ionized streamers and the molecular ring has been 
noted previously (Roberts & Goss 1993), and it is consistent with the suggestion that the western 
arc is the ionized inner edge of the molecular ring directly exposed to a central ionizing source 
(Telesco et al. 1996). Evidence for a similar link between neutral and ionized gas is also seen on 
the eastern rim of the CND, between the eastern arm and the clump I shown in Figure 5. From 
their comparison of the 8 GHz radio continuum map and the 30 ^m optical depth map, Telesco 
et al. (1996) postulated that the eastern arm is the ionized front of a protrusion of high-density 
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material facing the central ionizing source. The clump I identified in Figure 5 may indeed be this 

previously unseen high-density material, revealed for the first time. The extinction map derived 
from the HST/NICMOS Pa observations by Scovillc ct al. (2003) shows a hole in the obscuring 
material towards the northern arm and western arm, but high extinction regions with Ay > 30 
surround this hole, including the HCN clump corresponding to our clump I. Possible ionizing source 
candidates range from the putative central engine Sgr A* to hot UV stars near the Galactic center 
(sec discussion by Zhao et al. 1999). A more exotic possibility suggested by Macda ct al. (2002) 
is that the central source was much brighter in X-ray some 300 years ago as a result of the accretion 
of infalling material pushed in by the forward shocks from the supernova explosion in Sgr A East, 
leading to the ionization of the streamers. 

A close physical link between the neutral gas in the CND and the ionized gas in Sgr A West 
is a widely occurring phenomenon, and this fact can be demonstrated more convincingly in the 
comparison of the channel maps as shown in Figure 3. Within the velocity ranges of ±143 km 
s~^, the bulk of HAla and HCO"*" features trace the same rotation-like kinematics about Sgr A*. 
More importantly, the H41a features frequently occur just inside of the HCO"'" features (see the 
V = —91, —39, and -|-65 km channels). As in the comparison of the velocity integrated images 
(Figure 5), this physical link is more evident along the western arc. There is a tendency for the 
H41q; features to show an azimuthal offset from the HCO"*" features, particularly for the features 
closer to Sgr A* (e.g. V=— 39 and V=-|-65 km channels). These azimuthal offsets and a broad 
total line width (±200 km s~^) for the H41q; line are naturally explained if the gas follows nearly 
Keplerian orbits around a highly concentrated mass distribution, (see the discussion in § 4.5). 



4.5. Modeling Gas Kinematics in the CND 

Constrained by excellent high spatial resolution imaging data in 12.8 /ini [Nc II] and radio 
recombination lines, several plausible and detailed dynamical models for the ionized streamers 
have been constructed (e.g. Sanders 1998; Vollmer &: Duschl 2000). Models of gas kinematics 
for the CND have also been constructed, mostly to explain the features identified in low resolution 
images (^ ~ 10 — 20") that are susceptible to confusion with ubiquitous foreground and background 
emission (e.g. Liszt et al. 1985; Marshall et al. 1995). Taking advantage of the new high 
(and matched) resolution spectroscopic imaging data in H41q; and HCO"*", we can now examine in 
detail the gas kinematics and the gravitational potential in the central few parsec radius region in 
a self-consistent manner. Wc have already concluded that some of the ionized gas is orbiting Sgr 
A* in the plane defined by the CND (see § 4.4). Using the velocity information, we can also isolate 
and examine other kinematic features proposed to be associated with the CND, such as the "+50 
km s~^ streamer" and the "+70 km feature" (see Jackson et al. 1993). 

If a massive central source (2.5-4.1 xlO^ Mq, Genzel & Eckart 1999; Ghez et al. 2000, 
2003) dominates the gravitational potential in the central parsec region, the kinematics of the 
surrounding gas disk should be regular and well behaved. The neutral and ionized gas inside the 
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CND contributes only about 1000 Mq while we estimate the total mass of the CND to be about 
3 X 10^ Mq (see § 4.4). The mass of the central stellar cluster is more substantial and is highly 
centrally concentrated, doubling the enclosed total mass near 2 pc radius (Gcnzel et al. 2003). We 
adopt a mass distribution of the form similar to Eq. 3 of Vollmer &; Duschl (2002), after updating 
it using the most recent mass estimates of the massive central source by Ghez et al. (2003) and 
the dynamical mass estimates for the central 10 pc radius by Genzel et al. (2003), 

M{R) = (4.1 + 1.4i?^-25) X 10^ Mq (5) 

where R is radius in pc. The resulting Keplerian-like rotation of the gas disk nicely reproduces 
the observed features in the channel maps such as the azimuthal displacement for radially distinct 
features as discussed in § 4.4. To examine whether our interpretations of the gravitational potential, 
the geometry of the CND, and the relationship between the ionized streamers and the dense neutral 
gas in the CND are correct, we have created simulated channel maps of the CND using a Monte 
Carlo realization. They are shown side-by-side in comparison with the observed channel maps in 
Figure 11. 

Using Matlab software package, we created a disk of 3000 particles with a r~^/^ radial distri- 
bution. The model disk has a radial extent from 0.1 pc to 2.5 pc and has a radially dependent 
thickness with an opening angle of 4°. This particular disk geometry is adopted to allow a finite 
thickness which may play a role in the apparent gas distribution and kinematics in the projected 
view, but it is not motivated any particular observational constraints.^ The model disk is viewed 
with an inclination angle of 68° and is rotated by 60° in the sky plane to match the aspect ratio 
and the position angle of the CND. Each disk particle is in rotation about the central potential 
described by Eq. 5 and has a Gaussian random velocity component of 10 km s~^. To produce 
the model channel maps shown in Figure 11, model disk particles with the line-of-sight velocity 
matching the particular channel velocity window are selected, and their distribution is convolved 
with the observing beam (7.0" x 3.5", PA= —5°). The only adjustment made to the model was 
adding +13 km to the systemic velocity in order to match the highest velocity HCO"'" features 
at -117 km s~^ and +143 km s~^. 

Even though our model disk is constructed for an illustrative purpose only and is not tuned 
to match the observations in detail, the side-by-side comparisons of the model and data shown in 
Figure 11 offer significant insights on the gas distribution and kinematics: 

(i) As noted already, some of the ionized gas traced in H41a is sharing the same orbital plane as 
the dense neutral gas in the CND (see §4.4). The model disk has no knowledge of temperature 
or ionization information, but presumably the inner part (i? < 30") of the disk is either atomic 
or ionized. The clumps of the H41q; emission shown in contours on both panels in Figure 11 are 
frequently overlaying the model emission regions shown in grey-scale, especially in the channels with 



^The only estimate found in the literature is by Latvakoski et al. (1999) of ~ 7° opening angle, based on their 
own morphological and dust heating model. 
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velocities between it 91 km s~ . The northern arm is sometimes modeled as being displaced from 
the plane defined by the CND (e.g., Latvakoski et al. 1999), but a large fraction of the northern 
arm gas seen in the velocity range between and +143 km s~^ closely follow our kinematical 
model. 

(ii) The same comparison also shows that the gas in the eastern arm and the central bar do not 
follow the model. This leads us to believe that the ionized gas in the bar and the eastern arm is 
largely concentrated in a different orbital plane as suggested by others previously (e.g., Vollmer Sz 
Duschl 2000). 

(iii) The similarity between the observed HCO"*" emission and the model disk emission is more 
striking. Except for the channels with —39 km s~^ < V < +39 km s~^ where foreground absorption 
has severely altered the appearance of the emission features, most of the HCO''" emission seen in 
the channel maps have matching features in the model. The HCO"*" emission arising from the CND 
is frequently asymmetric about its major axis, unlike the model, and the CND is probably not 
a solid ring or a solid disk and is likely highly clumped. Two channels show substantial HCO"*" 
emission not predicted by our model, and these channels correspond to velocities for the "+50 
km streamer" and the "+70 km s~^ feature" that are already known to deviate from the CND 
rotation (see below). 

The +50 km s"^ streamer is a band of emission running nearly east- west at a slightly oblique 
angle in the V = +39 km s"^ and V = +65 km s"^ channels, previously noted by Jackson et al. 
(1993) and others. This feature extends well beyond the field of view of the OVRO interferometer 
and is several arcmin in length in the single dish HCN (3-2) channel maps of Marshall et al. (1995). 
This is clearly a structure much larger than the CND in size, and thus the +50 km s~^ streamer 
is most likely a foreground or a background structure, not directly associated with the CND. 

The +70 km s~^ streamer is the "7" shaped feature present in the V = +91 km s^^ channel. 
This structure is also seen in the single dish HCN (3-2) channel maps of Marshall et al. (1995), 
but it is more compact than the +50 km s"^ streamer. Although this may be a chance projection, 
this feature lies nearly exactly over the CND. However, its velocity differs from the CND rotation 
by about 130 km s~^. It is not as certain that the +70 km s~^ streamer is another foreground 
structure. Its long, linear appearance is suggestive of a large scale shock in origin. It is not seen in 
H41a or in continuum to offer any further clues, and the lack of any ionized gas at the corresponding 
velocity is an important argument against the +70 km streamer being part of the CND. 

One of the key arguments against the CND being a coherent disk or a ring is the apparent 
gap along its eastern side in the previous HCN and HCO+ maps. This gap now appears filled 
with several compact gas clumps such as the clump I in our new HCO+ image shown in Figure 5. 
In addition, neutral gas with higher temperature and density also fills this gap as inferred by the 
observations of [C II] and NH3 (see Jackson et al. 1993; Herrnstein & Ho 2002). Therefore earlier 
arguments based on the apparent morphology in the lower sensitivity HCN and HCO+ maps are 
no longer valid. The mean velocity of the warm molecular gas found by Herrnstein & Ho (their 
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clump D) is slightly redshifted (< V >~ +50 km s~^) than expected from our simple rotation 
model, but this velocity shift is in the correct sense if their average spectrum includes the emission 
from the entire north-south extent of their clump D. 

In summary, we have successfully demonstrated the continuity in the distribution and kine- 
matics of the ionized and the neutral gas within the central 3 pc diameter region surrounding 
Sgr A* using our new observational data and a simple kinematical model. Our kinematical model 
reproduces many of the observed neutral and ionized gas features, as shown in Figure 11. The 
success of this simple model strengthens the notion that most of the gas found in the central parsec 
of Sgr A* partake in a simple rotation around the central source within the disk plane defined by 
the CND. Our kinematical model also identifies very clearly which gas components do not follow 
this simple rotation. In addition, our model suggests that some of the subtle observed effects such 
as the slight asymmetry between the near and far side of the disk are the products of the finite size, 
asymmetric observing beam and projected viewing of the 3-D CND structure (e.g., see V = ±39 
km s~^ channel). The overall picture of the gas distribution and kinematics in the central 3 pc 
surrounding Sgr A* is not yet complete because information on the distribution and kinematics 
of atomic gas inside the CND is quite incomplete. Future arcsecond resolution imaging of dense 
atomic gas tracers such as [C II] will provide important clues for completing our understanding of 
the gaseous environment in this region. 



5. SUMMARY 

The OVRO millimeter array observations of the H41a and 92 GHz continuum emission from 
the ionized streamers around Sgr A* (Sgr A West) and the HCO"*" (1-0) emission from the circum- 
nuclear disk (CND) are presented along with a new survey of the hydrogen recombination lines at 
millimeter wavelengths (H31a, II41a, and II44a) obtained using the NRAO 12-ni telescope. These 
observations are compared with the VLA II92a and the 8.3 GHz continuum data by Roberts &; 
Goss (1993). We summarize our results as following: 

1. The 92 GHz continuum image is essentially identical to the VLA continuum image at 8.3 GHz 
in morphology. The total continuum flux detected inside the 85" diameter is about 12 Jy, but some 
extended flux associated with structure larger than about 30" is missed. The brightest feature, Sgr 
A*, is unresolved at 3" resolution of our observations and has a peak flux density of 1.95 ± 0.30 Jy 
beam~^. Assuming the 7000 K electron gas fills our 0.005 pc (1000 AU) beam, we infer average 
continuum optical depth of r < 10~^ and electron density rig ~ 2 x 10^ cm~^. 

2. The kinematics and the morphology of H41a emission are very similar to those observed in 
H92a and 12.8 fim [Ne II] lines. The brightest emission occurs along the north arm with a peak 
line brightness temperature of less than 10 K. This suggests a H41a optical depth of tl < 10~^ 
averaged over the 13 km s~^ channel width. The emission along the central bar displays a line- 
of-sight velocity ranging over 300 km s~^, and this gas may be located in the deepest part of the 
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gravitational potential within the CND region. 

3. The molecular CND is imaged in HCO"*" with a 40% improvement in resolution and 4-15 times 
better sensitivity over the previous BIMA observations made in HCN and HCO+. The CND 
consists of several discrete concentrations, and their peak brightness ranging between 10 and 30 
K averaged over the 0.1 pc synthesized beam suggests that they originate from the warm (60-100 
K) gas responsible for NH3 and dust continuum emission in the CND. A virial analysis suggests 
that these clouds have mean density of n ^ 10^ cm ^, sufficient to withstand the tidal shear in 
the region, and the CND may be a long lived structure. We estimate a total gas mass of about 
3 X IO^Mq for the CND. Our non-detection of HCO+ absorption along the Sgr A West places an 
upper limit to the gas column density of < 4 x 10^^ cm~^, or the neutral gas within the CND is 
depleted of molecules. 

4. A line-to-continuum ratio analysis under the LTE assumption yields a fairly uniform electron 
temperature of about 7000 K along the ionized streamers in the Sgr A West complex, in an excel- 
lent agreement with the electron temperature calculated at 8.3 GHz using the H92q; observations 
(Roberts Sz Goss 1993). Evidence for increased (or an unusually low L/C ratio) along the 
central bar just southwest of Sgr A* is found in our data as well even though our data include all 
line emission within ±350 km s~^. The low L/C ratio we find at 92 GHz challenges the optically 
thick continuum sources IRS 2 and IRS 13 as additional sources of continuum emission as suggested 
by Roberts, Yusef-Zadeh, &: Goss (1996). Therefore, the elevated Tg (a lower L/C ratio) seen near 
the location of IRS 2 and IRS 13 is still not fully understood. 

5. A comparison of the velocity integrated line flux densities for H31a, H41a, H44q!, and H92a 
emission from the Sgr A West complex shows line ratios characteristic of spontaneous emission from 
optically thin gas near LTE. A detailed analysis reveals that observed recombination line brightness 
is reduced by 10-30% (i.e., 6„ ~ 0.7 - 0.9) from the LTE value. 

6. A combined analysis of the new high resolution H41q; and HCO"*" images and a kinematic model 
clearly demonstrates a widely spread, close physical and dynamical link between the molecular gas 
in the CND and the ionized gas, including along the eastern rim where a gap in the CND was 
previously reported. The same analysis also clearly identifies several features that do not follow the 
disk rotation. We conclude that the "+50 km s~^ streamer" is not directly associated with the 
CND. The morphology of "+70 km streamer" is suggestive of its association with the CND, 
but its kinematics and absence of any ionized gas challenge the idea of their close association. 

We are grateful to D. Roberts for providing us with the H92a and 8.4 GHz continuum data from 
the VLA for comparison. We are indebted to C. Lang for her expert advice and help with the GIPSY 
data reduction package and for many helpful suggestions as the referee of this article. In addition, 
we would like to acknowledge the direction gained from insightful discussions with W. M. Goss and 
late K. R. Anantharamaiah. We would like to thank the OVRO and the NRAO 12-m Telescope 
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Fig. 1. — Deconvolved image of 92 GHz continuum emission from Sgr A complex at 6.95" x 3.47" 
(PA = —5°) resolution of the synthesized beam. The cross marks the position of Sgr A*. The 
grey-scale range is between 0.03 to 0.40 Jy beam~^ with RMS noise of 0.01 Jy beam~^. 
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Fig. 2. — Integrated line emission of the H Ala hne from Sgr A complex at 6.95" x 3.47" (PA = 
—5°) resolution of the synthesized beam. The cross marks the position of Sgr A*. The grey-scale 
range is between and 12 Jy km beam~^. 



Fig. 3. — (a) Channel maps (from —221 km to —13 km s~^) of HCO"'" emission are shown 
in grey-scale and H41a in contours, (b) Channel maps (from +13 km s"^ to +221 km s~^) of 
HCO"*" emission are shown in grey-scale and H41a in contours. The grey-scale ranges between —100 

(white) and +500 (black) mJy beam~^ for HCO"*". The contour levels are —2, +2, +3, +4, +5, 
and +6 times 15 mJy beam~^ (Icr)- The location of Sgr A* is marked at the center with a cross. 
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Fig. 4. — A contour map of H41a emission is shown with hne-to-continuum (L/C) plots of respective 
regions. The contour levels are +1, +2, -|-3, +4, -|-5, and +6 times 2.3 Jy km s^^. Ratios of the 
L/C are marked on the y-axis and the velocity ranges in km s"-*^ are shown on the x-axis. Crosses 
mark the data points, solid line represents the model Gaussian fit, and the dashed line represents 
the residual. 
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Fig. 5. — Integrated line emission of the HCO'*' line from Sgr A complex at 6.95" x 3.47" (PA = 
—5°) resolution. The cross marks the position of Sgr A*. The grey-scale range is between 7 to 32 
Jy km s~^ beam~^. Superimposed are the contours of H 41a emission. The contour levels are 2.3, 
4.6, 6.9, 9.2, 11.6, and 13.9 Jy km beam~^. The boxes mark the dense gas clumps described in 
Table 5. 
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Fig. 6. — HSla, H41a, and H44a spectra of Sgr A West obtained with the NRAO 12-m telescope. 
The H31a spectrum is made by summing 8 independent spectra covering the whole Sgr A West, 
scaled down by a factor of four. Both the OVRO (solid line) and the NRAO 12-m (dotted line) 
spectrum are shown for H41q; line. The excess emission in the NRAO 12-m spectrum near -1-200 
km s~^ may be due the contamination by widely distributed foreground CH3CN emission (see text). 
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Fig. 7. — Grey-scale representation of LTE electron temperature distribution in Sgr A West. The 
grey-scale ranges between 3000 to 15000 K. The 8.3 GHz continuum image at 0.62" x 0.51" resolution 
by Roberts k, Goss (1993) is shown in contours. The contour levels are 2.5, 5, 10, 20, 30, 40, 50, 
and 60 mJy beam~^. The thick cross marks the position of Sgr A*. The location of the IR source 
IRS 2 is shown with a thick white arrow, and the second continuum source about 2" to the north 
is IRS 13. 



-29- 




Fig. 8. — A ratio map between H92a and H41q; line emission is shown in grey-scale (from 0.0 to 
1.5). The integrated H92q; line data from Roberts et al. (1993) were convolved with the OVRO 
beam. The ratio is nearly constant, ~0.15, in most regions. 
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Fig. 9. — Velocity integrated line flux densities of H31a, H44a, and H92q! are plotted as a function 
of frequency. The H41q: measurement from the OVRO data is shown only as an upper limit since 
some of the extended emission is resolved out by the interferometer. The solid line is the best fit 
power-law relation, Sl oc v°'. The slope a = 0.95 ± 0.18 is very close to the theoretical value for 
spontaneous emission {Sl oc 
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Fig. 10. — A comparison of the H41a (solid line) and H92q! (dotted line, multiplied by a factor of 
15) spectrum, one at the Sgr A* position (right) and at an adjacent position 7" to the east (left). 



Fig. 11. — (a) Channel maps (from —195 km s~^ to —13 km s~^) of HCO^ emission are shown in 

grcy-scalc. In each twin panel plot, the observed HCO^ emission is shown on the left panel while 
the simulated disk model (see § 4.5) is shown on the right panel. The contours of H41a emission are 
shown in both panels. The grey-scale ranges between —100 (white) and +500 (black) mJy beam^^ 
for HCO"*". The greyscale for the model shown on the right panel uses the full range in model units. 
The contour levels are —2, +2, +3, +4, +5, and +6 times 15 mJy beam~^ (lo")- Each box is 100" 
on each side, and Sgr A* is located at the center, (b) Channel maps (from +13 km s~^ to +195 
km s~^) of HCO+ in grey-scale superimposed with contours of H41a. 
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Table 1. Summary of the OVRO Observations 



Rest Frequency of H41a line 
Rest Frequency of HCO+(J = 
Primary beam (HPBW) 
Phase center 

Center LSR velocity 
Synthesized beam (HPBW) 
Spectrometer settings 
RMS noise 



92.034 GHz 
1 0) line 89.188 GHz 

85" 

q;(51950) = 17'*51"*29.4* 
5(51950) = -28°59'19.2" 

0.0 km s"^ 
6.95" X 3.47" (PA= -5°) 
60 X 4 MHz (13 km s-^) 
0.13 Jy beam-i (0.21 K) 



Table 2. 92 GHz continuum emission 



Feature 


Size 






Stotal 




(pc) 


(Jy Beam~^) 


(K) 


(Jy) 


Northern Arm 


0.9 


0.39 ±0.07 


3.1 ±0.6 


1.09 ±0.16 


Sgr A* 


0.1 


1.95 ±0.29 


15.6 ±2.3 


1.95 ±0.30 


Bar 


1.1 


0.75 ±0.14 


6.0 ± 1.1 


0.90 ± 0.29 


Western Arc 


1.5 


0.14 ± 0.02 


1.1 ±0.2 


1.08 ± 0.08 


Eastern Arm 


0.4 


0.29 ± 0.05 


2.3 ±0.4 


0.52 ±0.12 


Table 3. 


A S 


ummary of H41 


a line intensity 


Region 


Size 


Speak 


Tpeak 


Stotal 




(pc) 


(Jy Beam~^) 


(K) 


(Jy) 


Northern Arm 


0.9 


1.18 ±0.27 


9.4 ±2.2 


4.70 ±0.71 


Bar 


1.1 


1.12 ±0.21 


9.0 ± 1.7 


5.70 ±0.37 


Western Arc 


1.5 


0.50 ±0.14 


4.0 ± 1.1 


1.86 ±0.39 


Eastern Arm 


0.4 


0.90 ±0.15 


7.2 ± 1.2 


2.02 ± 0.36 
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Table 4. LTE electron temperature of the regions identified in Figure 4 



Region 


Box Center 


Vlsr 


Ay 


< Tl/Tc > 


Temperature T* 




(J2000) 


(km s~^) 


(km s~^) 




(K) 


1 


17'^42"^29.^71, - 


28°59'09."2 


+114 ±3 


78 + 5 


0.50 + 0.04 


6394+608 


2 


17''42"^28.-'60, - 


28°59'18."2 


-132 ±8 


213 + 18 


0.18 + 0.02 


6473+891 


3 


17''42"^27/'92, - 


28°59'46."7 


-92 ±3 


43 + 5 


0.80 + 0.13 


7470+1372 


4 


17^42™29.'^27, - 


28°59'19."2 


-117 ±20 


347 + 48 


0.03 + 0.01 


18806+4770 


5 


17''42"^30/32, - 


28°59'26."2 


+142 ± 3 


70 + 5 


0.66 + 0.07 


5528+623 


6 


17'^42"*29.^59, - 


28°59'17."2 


+46 + 3 


101 + 6 


0.26 + 0.02 


9057+825 



Table 5. HCO"*" line intensity of the nine brightest clumps in the CND 



Clump 


Size 


Speak 












(pc) 


(Jy Beam~^) 


(K) 


( km s~^) 


(Mo) 


(cm~^) 


A 


0.4 


3.5 + 0.8 


27.9 + 6.4 


20.7 


2.1 xlO^ 


1.4 x 10^ 


B 


0.3 


1.6 + 0.4 


12.7 + 3.2 


17.3 


1.1 xlO^ 


1.8 X lO'^ 


C 


0.3 


1.0 + 0.8 


8.0 + 6.4 


12.0 


0.5 xlO"^ 


0.9 X lO'^ 


D 


0.4 


1.2 + 0.3 


9.5 + 2.3 


29.9 


4.5 xlO"^ 


3.0 X lO'^ 


E 


0.2 


1.5 + 0.9 


11.9 + 7.2 


11.6 


0.3 xlO^ 


1.8 X 10'^ 


F 


0.3 


1.9 + 0.4 


15.1 + 3.2 


19.4 


1.4 xlO^ 


2.2 X 10'^ 


G 


0.6 


1.5 + 0.4 


12.0 + 3.2 


16.1 


1.9 xlO'' 


0.4 X 10^ 


H 


0.4 


1.0 + 0.5 


8.0 + 4.0 


19.0 


1.8 xlO^ 


1.2 X lO'^ 


I 


0.2 


1.8 + 0.7 


14.3 + 5.6 


13.0 


0.4 xlO"^ 


2.3 X lO'^ 



-34- 



Table 6. Summary of the NRAO 12-m Observations 



Line 




SAV 




(GHz) 


(Jy km s^-"^) 


H31a 


210.502 


1034 ± 207 


H35a 


146.969 




mia 


92.034 






74.645 


232 ± 46 


H92a 


8.309 


42.6 ±6.5 
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